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In optimizing the topology of wireless networks built of a dynamic set of spatially embedded agents, there ai'e 
many trade-offs to be dealt with. The network should preferably be as small (in the sense that the average, or 
maximal, pathlength is short) as possible, it should be robust to failures, not consume too much power, and so 
on. In this paper, we investigate simple models of how agents can choose their neighbors in such an environment. 
In our model of attachment, we can tune from one situation where agents prefer to attach to others in closest 
proximity, to a situation where distance is ignored (and thus attachments can be made to agents further away). 
We evaluate this scenario with several performance measures and find that the optimal topologies, for most of the 
quantities, is obtained for strategies resulting in a mix of most local and a few random connections. 



I. INTRODUCTION 

The performance of a network is a consequence of three 
factors — the hardware, the network protocols and the net- 
work topology. For some networked systems, like the wire- 
less mesh networks, the topology can be controlled rather 
easily through the medium access layer, and it is thus practi- 
cally possible to optimize the topology of the network. In this 
work, we investigate a scenario of static agents (sensors, indi- 
viduals with Wi-Fi devices, etc.), embedded in space, joining 
and leaving the system (so called "churn"). Our scenario is 
based on the following three assumptions. First, only new 
agents and neighbors of an agent leaving the network cre- 
ate new links. Second, the power consumption of each agent 
is, for fairness and practical reasons, to be restricted. Third, 
the agents are localized with uniform randomness over the 
unit square. From these starting points, we evaluate differ- 
ent simple attachment strategies to optimize several network- 
structural quantities, capturing efficiency of communication, 
scalability and robustness to failure. 

Many of the previous studies on topological optimization of 
wireless networks have concerned static point patterns (TT), 
not (as in our study) systems under churn. Nevertheless, these 
studies form a theoretical backdrop of ours. The simplest 
approach to connecting wireless agents in the plane is to let 
all agents within range be connected. Assuming identical 
agents, one arrives at the Unit Disc Graph model of wire- 
less networks (19). Maintaining and routing information in 
a network with this architecture, however, is a waste of power 
and memory resources. Many theoretical efforts in this area 
have been focusing on energy efficient, yet scalable topolo- 
gies — a seminal work is (l23l l. where the "Yao graph" was 
proposed. This graph is constructed by connecting each node 
to its closest neighbor in k equiangular sectors. It has sev- 
eral topological features, such as (almost surely) full connec- 
tivity and scalable power consumption, which are desirable 
for wireless networks. Subsequent studies have proposed im- 
proved proximity graphs, with localized construction algo- 
rithms more til3\ [TH or less (|9j [TOl [TSj) similar to the 



Yao graph, and studied general features of minimum power 
topologies (il8i i21}. Our approach is from a slightly differ- 
ent angle, taking more inspiration from the recent complex 
network literature ([3] IS [171 1. These works discuss how dy- 
namic properties of networks in general (including commu- 
nication performance in technological networks) are related 
to network structure (how the observed network topology dif- 
fers from what is expected in a unconstrained random network 
model). This approach lends itself naturally to analyzing top- 
ics like point-to-point communication, and data distribution 
scenarios in large, to some extent random, systems. In this 
work we evaluate the possibility to optimize network topology 
in the case when connections cost power, and consequently, 
due to limitations of energy storage, agents maintaining long 
connections can keep only a limited number of neighbors. 

In the next section, we give a precise description of the 
problem and our proposed attachment models and discuss re- 
lated practical issues. We evaluate the proposed models con- 
sidering a large set of performance measures in Section III 



Finally, Section IV concludes our work. 



II. PRELIMINARIES 

A. Problem statement 

We consider a set V of agents spread out on the unit 
square. In our numerical study, the spatial distribution will 
be uniformly random, but generalizing to other point patterns 
is straightforward. At every time step, one random agent 
leaves the system and a new one enters (at a random coor- 
dinate). The agents form a graph with bi-directional con- 
nections. For two agents / and j at coordinates r, and to 
be connected it costs both agents a power consumption of 
p{i, j) - |r,- - Tjf {2 < 6 < 4) dSj. In this work we focus 
on the limit 6-2. We assume the agents to be selfish and ho- 
mogeneous in the sense that no agent / accept a higher power- 
consumption Ejsr, Pii' j) (where F,- is /'s neighborhood), than 
^'max- This also ensures scalable, total power consumption. 
We also assume a lower power limit Pmin, such that while an 



2 



agent consumes less power than f min, it will attempt to create 
new links in the network (unless creating any new link would 
result in a larger power consumption than Pmax)- Thus P^i„ 
represents a minimal investment all agents agree to contribute 
with to the system. Based on these precepts, we try to find 
simple attachment strategies (rules for agents to select other 
agents to link to) that optimize various objective functions. 
To assess unperturbed network performance, we measure the 
following quantities (that will be discussed in greater detail in 
the Section|lIl|. 

• Connectivity. We measure connectivity as the fraction 
of time steps when the graph is completely connected. 
This quantity should, in any functioning network, be 
very close to one. Except when we explicitly study this 
quantity, we tune the parameters so that the graph is 
almost always connected. 

• Power consumption and power efficiency. We con- 
sider the value of aggregate power consumption and the 
power efficiency (the ratio of power consumption across 
the shortest multihop path and a direct link). Low val- 
ues are desirable for efficient communication. 



B. Model definition 

We will investigate two models,' both having in common 
that they interpolate between situations where the power bud- 
get is spent only on short-range connections, to a situation 
where connections are made at random regardless of location. 
At every time step, a random agent /' is removed along with 
its connections F,' and a new agent added. This can cause 
both the new agent and F,' to have power consumptions below 
^'min- All agents / with P(i) < Pmin try to attach to new agents 
j. If adding a link (/, j) would make either P{i) or P(j) larger 
than PmsLx, or if (/, j) akeady is a link in the graph, the link 
is not added. These properties are the same for both models. 
The differences between the models are: 

• Model A. With probability q^, attach / to the closest pos- 
sible agent, otherwise attach / to a random agent. 

• Model B. There are two kinds of agents — a fraction qB 
of them always attach to the closest possible agent, the 
rest attaches to a random agent. 

Note that in the limits qA - qs - ^ and q^ - qg = 1, the 
models are identical. 



Diameter dmax - max/j-gv d(i, j). Where d{i, j) is the 
distance (the number of links in the shortest path) be- 
tween / and j. 



Average distance d 



2; jev d(i, j). The over- 



head in communication time grows with distance; thus, 
a small diameter and short average distance are desir- 
able. 



Spectral gap, E — the difference between the largest 
and second largest eigenvector of the adjacency matrix. 
This is an approximate measure of the link expander 
property (l8l) that embodies several desirable properties 
for wireless networks (little redundancy, high robust- 
ness, efficient broadcast, etc.). 



Since we do not construct an optimal network each time 
step, only update the network incrementally, it is important to 
monitor the effect of churn. Indeed, this problem is similar to 
another important aspect — the robustness to failure and we 
measure it by the two quantities: 

• The change, Ac/^ax. of the diameter as a random node is 
deleted from the network. 



• The change. Ad, of the average distance. 

We aim at finding simple, local attachment strategies opti- 
mizing these quantities — local in the sense that each agent 
follows its own rules. 



C. Practical considerations 

The main assumption of the proposed models is that 
the power consumption of an agent / is proportional to 
Yjjer, P('' y)' where F, is the set of agents connected to / and 



p{i, j) 



This requires transmission power control 



to be implemented at the agents. We motivate the assump- 
tion on power consumption as follows. First, once the net- 
work topology is defined, and if the traffic load in the network 
is balanced, both unicast and broadcast routing scenarios are 
likely to have a similar average usage of the connections. (In 
a situation with heterogeneous agents of very different capac- 
ities and network positions, this would not hold.) Second, 
considering dynamic power control, the required transmission 
power depends on the transmission channel characteristics. In 
systems designed for tolerable outage probability it is suffi- 
cient to consider the distance-dependent path loss to set the 
required transmission power ( 18). Finally, we do not consider 
the energy spent for data reception and processing. These are 
distance independent additive values and consequently would 
not change the main conclusions of our study. 

For the implementation of the attachment strategies de- 



scribed in Section II. B the problems of distance-dependent 



neighbor selection and efficient access control to the shared 
wireless medium have to be solved. 

If the network, for efficient higher layer protocols (like ge- 
ographic routing (20)), keeps information about the location 



A Java Applet implementation of the model can be found at http://| 



statphys . skku . ac . kr/~bjkim/ Applet/mobile .html 
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FIG. 1 The average degree k for our models, (a) shows results for 
model A, (b) for model B. = '^Pmin = 2, and (c) shows the 
size scaling for qA = 1b = 0. The line in (c) represents growth 
proportional to W''-. Standard errors are smaller than the symbol 
size. 



of the nodes, then such information can also improve the pro- 
posed attachment models. If absolute node location informa- 
tion is not maintained, agents can select neighbors by discov- 
ering a given area via beaconing with increasing transmission 
power (see for example (21) for detailed solutions). 

The medium access control is responsible for sharing the 
wireless medium in an efficient way. TDMA, CSMA and 
mixed solutions are usually proposed, depending on the ex- 
pected network load (16). Since the interference region of a 
connection depends on the distance of the transmitting and re- 
ceiving agent, the construction of an optimal access control 
scheme is not trivial in the addressed scenario. An optimal 
access control should take the attachment model into account 
and schedule most local and random connections separately. 
The efficiency of the medium access control can potentially 
affect the network performance, which calls for extended fu- 
ture studies. 



III. NUMERICAL RESULTS 

In this section we will discuss results from our numerical 
simulations. Unless otherwise stated, we use 10'' time steps 
for equilibration and 10^ time steps for measuring averages. 



A. Average degree 

Many of the performance measures, like connectivity and 
network diameter are directly related to the number of con- 
nections the agents maintain — that is, the average degree. 



Since the number of connections an agent can maintain in the 
considered case depends on the attachment strategy, we start 
investigating the average degree. 

We can derive approximate results for the two extreme 
cases, qA - iB - and - qs - ^- For qA - qs - 0, each 
agent tries to connect to the closest available neighbor. As- 
suming that all connections can be made and that the agents 
are uniformly spread out (so that the number of agents 
within distance r from point is nr^N), this means that the 
power consumption to attach to agents within a distance R is 



2nN 

R 

2 + 6 



2+i 



(1) 



For a given minimal power Pmin the number of agents within 
reach is 



^9 + A ^ 2/(2+5) 



(2) 



or in other words, for 6-2 the average degree k e 0(N^^^). 
This relationship is confirmed in Fig.[T] In the opposite limit, 
qA - qs - 1^ all agents attach to others at random. Since the 
average distance to another, random, agent is independent of 
A^, then so is k (which also can be seen in Fig.[T]l. The func- 
tional forms of k, as functions of qA,B are strikingly different. 
Even if k itself is not a performance measure, it affects net- 
work performance; therefore a stable response to the param- 
eter values is beneficial. We note that the curves for Model 
B has a smaller slope for small values of qs, whereas Model 
A has larger magnitude of the derivative when qA is close to 
zero. 

The key to understanding the difference in degree between 
the two model is the observation that the number of long-range 
links is much lower in Model B than in Model A (Fig.[T]). P^^^ 
controls the power consumption so that it is relatively similar 
for both models. With the power consumption constrained, 
more power-expensive, random links leads to a sparser net- 
work. So why does Model A have more random links? New 
links are added until the power consumption reaches above 
Pmin- Since the steps in power consumption are larger when 
a random connection is added, the chance that the last added 
link, for Model A, is by random attachment is larger than qA- 
For model B, on the other hand, the chance the last added con- 
nection is long range is qg, so k will be a linear combination 
of the qs - and q^ = 1 limits. Therefore, for qA - qs the 
number of long-range edges and the power consumption will 
be higher for Model A than Model B. These mechanisms can 
also be understood from our example figure. Fig |2] To make 
this figure readable we use smaller values of and P„in than 
in the rest of the paper. We also use a relatively large Pmax, 
making the k more similar between the models at the expense 
of a more different power consumption. From this figure we 
see the difference in number of long-range links even more 
pronounced than in Fig.[T] 

From Fig.[T|a) and (b) we learn that the degree can be con- 
trolled continuously by the model parameters qA and qs- It 
would maybe be convenient for the discussion if we could 
neglect the model parameters and investigate the models in 
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FIG. 2 Snapshot of the models with N = 100, Pmin = 0.1, A^ax = 1, 
and = qs = 0.1. The dashed lines are formed by random attach- 
ment. The triangles in (b) are the agents with a random attachment 
strategy rather than attachment to the closest other node. 



FIG. 4 Power consumption for the two attachment models A and B 
as functions of their respective parameters, q^ and q^. Other param- 
eter values are the same as in Fig.[T] Standard errors are smaller than 
the symbol size. 
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FIG. 3 The connectivity is shown as a function of q^ g and P,nin 
(panel (a) and (b)) or P,nax for Model A (panel (a) and (c)) and Model 
B (panel (b) and (d)). The fraction ip of time steps when the network 
is fully connected is measured and displayed as color scales so that 
the numbers at the color bars correspond to - lg(l - (p), white means 
that the network is almost always connected. In all panels we set 
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and use N = 1000. 



B. Connectivity 

A fundamental functional requirement is that the network 
should be connected. There is, in principle, nothing that guar- 
antees connectivity in our construction algorithm. We inves- 
tigate the connectivity, quantified as (f — the fraction of time 
steps when the network is fully connected. Since (p is close to 
one, we rather plot - lg(l (Fig.[3]l. The larger are values of 
the energy limits (both Pmin and Pmax), the better is the con- 
nectivity — by consuming more power, the network can al- 
ways be made connected — a trivial but necessary fact. In the 
white region of the density plot in Fig. [3j the network is fully 
connected more than 99.9% of simulation runs {(p > 0.999 
and thus - lg(l - 1^) > 3). In the rest of this paper we will (un- 
less stated otherwise) investigate a region of parameter space 
where the network is almost surely connected — f = 1 and 
P -7 

^ max — 

As seen in Fig. [3] a qualitative difference between the two 
models is that the connectivity of Model B is less dependent 
on qs than Model A is dependent on q/^. For a large range 
of qs (Fig. [3];b)), the network is connected almost regardless 
of Pmin- One explanation is the higher degrees of Model B 
— adding links, increasing the degree, can connect, but not 
fragment, a network. The connectivity depending on the P^ax 
value (Fig. [3|c) and (d)) shows different tendencies: Model 
A becomes connected at lower f ^ax-values (and consequently 
lower Pmin-values). To explain this feature we evaluate the ac- 
tual power consumption of the two models in the next section. 



terms of some basic network quantity like k. However, such 
an approach is complicated by the fact that, for ^a,b < 1, the 
average degree also depends on (see Fig.[TJc) and the above 
discussion). Moreover, the size scaling of k depends on qA and 
qs, going continuously (but with different functional forms) 
from k ~ N^l^ for qA - qs - 0, to size independence for 
qA - qs - ^- For this reason we keep the control parameters 
qA,B as our main independent parameters in this investigation 
and do not go into details of functional forms when discussing 
the size dependence of performance measures. 



C. Power consumption 

The power consumption of an agent is a stochastic variable 
in the interval [0, Pmax], changing as agents join and leave 
the network. (While in both models the agents try to allo- 
cate at least f min, we have to note that if an agent cannot 
connect without raising the power consumption over Pmax, 
its power consumption will be below Pmin-) If the average 
power consumption is higher, so more resources are invested 
in the infrastructure, one can expect the network properties to 
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FIG. 5 The average (per node pair) ratio between tlie actual con- 
sumed power and the power of a direct connection p for the two at- 
tachment schemes as a function of their respective control parameter 
Other parameter values are the same as in Fig.[T| Standard errors are 
smaller than the symbol size. 



improve. Therefore, to better be able to analyze our perfor- 
mance measures, we investigate the level of power consump- 
tion. Fig.|4|a) and (b) shows the average power consumption 
with f = 1 and f ^ax - 2. The power consumption max- 
imally fluctuates ~ 20% and it is always lower for model B. 
This difference increases with system size for all values ex- 
cept close to qA,B = 1, meaning that performance advantages 
for Model A in these limits, especially for larger systems, 
should be evaluated with this increased power consumption 
in mind. Or, alternatively, a specific g'^i-value can be trans- 
lated to a (higher) qB-valne, corresponding to the same power 
consumption. 

Why is flie power consumption higher in Model A? Again, 
we use Fig. |2] can give an illustration. Above, we explained 
why, if qA = qs. Model A has larger fraction of long-range 
edges than Model B. Since the steps in power consumption is 
larger when adding random edges, if a random edge is the last 
one added, then the power consumption overshoots Pmi,, more. 
In sum, more randomly added links mean more overshoot and 
thus higher power consumption. 

Power efficiency in multihop networks is often measured 
by the ratio p of actual consumed power to the power needed 
for a direct connection, averaged over all vertex pairs (scaling 
like the distance squared): 



- (N\ Zj 

[2) i<jeV 



P(i, j) 



(3) 



where n{i, j) is the set of links forming the shortest path (in 
number of hops) between / and j (19\. If the route between 
the source and the target is circuitous, p might be much larger 
than one, but if it proceeds in several steps relatively straight 
towards the target, p can be smaller than unity. Indeed, the 
lower bound of p is zero — assume / and j, at distance r 
are connected by a straight path of n equidistant links. Then 
/?(/', /) for all links along the path is r^/ir^, so the whole sum 
in the numerator equals r^/n, while the denominator is r^, giv- 
ing a contribution 1 /« from the node pair {/, j} to p. In Fig. |5] 
we plot this quantity for Model A and B, and two different 
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FIG. 6 Distance-related measures for different sizes and parameter 
qA,B values. Panels (a) and (b) shows the average distance d while 
(c) and (d) shows the diameter (i„,ax- Other parameter values are the 
same as in Fig.[T] Standard errors are smaller than the symbol size. 



network sizes. The p-curves are monotonically increasing for 
both sizes and attachment models. Both models have some 
regions of their parameter space where p < 1, which is the 
definition for being a power spanner (W). Model B shows 
lower values of p for a large region of its parameter space, 
which means that shortest path routes are straighter compared 
to the ones in Model A, a results of the different distribution 
of long-range random connections. 



D. Distances 

Both for unicast and multicast communication, short dis- 
tances (in number of hops) are beneficial, as they lead to lower 
delays (4). In Fig.|6]we investigate the average distances and 
diameter, corresponding to mean and extreme-case connec- 
tions as a function of the model parameters qA and qs- For 
both Model A and B, and both distance measures, the minima 
are attained for intermediate values of qA^s- In other words, a 
mix between random and short-range connections optimizes 
network distances. Furthermore, the minimum occurs for low 
?A,B-values, i.e. with larger proportion short-range connec- 
tions compared with the random connections. The decreasing 
distances for qA.B ~ can be understood as a "small-world ef- 
fect" — we only need to introduce a few random links to a reg- 
ular graph for it to go from algebraic to logarithmic distance 
scaling (fni22l). The increasing distances for larger ^/i,B-values 
can be explained by the decreasing average degrees. Compar- 
ing Models A and B, we note that Model A gives the smallest 
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FIG. 7 The spectral gap, the difference between the two largest 
eigenvalues of the adjacency matrix, for model (a) and (b). Other 
parameter values are the same as in Fig.[T] Standard errors are about 
the size of the symbols, but omitted for readability. 
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FIG. 8 The expected difference in average distance as a random node 
is deleted from the network. Panels (a) and (b) shows curves for 
Model A and Model B respectively. Other parameter values are the 
same as in Fig.^ Standard errors are smaller than the symbol size. 



distances (both average and maximal). If one increases the 
degree of a graph by adding links to it, its average distance 
cannot increase. Similarly, most network models have dis- 
tances decreasing with the average degree. With that in mind, 
it is a little surprising that Model A has both shortest dis- 
tances and smallest average degree. This non-trivial, purely 
topological effect can only be related to the distribution of 
longer-range, random links. In Model A, the long-range links 
are more evenly distributed out between the agents (compared 
with Model B where, at least one side of the every long-range 
link belong to one class of nodes). Model A loses this advan- 
tage as more random links are inserted in the network, then 
the average and maximal distances increase together with the 
decreasing degree according to Fig. [T] Model B has the ad- 
vantage that it is not as sensitive to qs as model A is to qa- 

E. Expander property 

For multicast and broadcast communication, redundancy is 
an important topological factor. Imagine a triangle between 
agents 1, 2 and 3. A broadcast originating from 1 reaches 2 
and 3 in two steps. Then, however, the link between 2 and 
3 is redundant — both 2 and 3 already have the packet. If a 
graph has few such redundancies, it is said to have good ex- 
pander properties. The usual way of defining this property is 
via the expansion factor <t) — the minimum over all subsets 
S , smaller than half the entire vertex set, of the ratio between 
the size of the neighborhood of S and the size of S (5). <1) is a 
computationally hard quantity. Instead, we measure the spec- 
tral gap E — the difference between the leading and second 
largest eigenvalues of the adjacency matrix, which is known 
to be a lower bound to O (apart from a factor 2) 

In addition to efficient broadcast communication, large 
spectral gap also leads to two other desirable network prop- 
erties. First, it reflects robustness. It is known that if the graph 
has a large spectral gap one needs to cut many links to split it 
into two disjoint subnetworks (8). Consequently, the spectral 
gap is a measure of resilience against a worst-case scenario of 
an adversary deliberately trying to disconnect different parts 



of the network 0. Second, networks with large spectral gaps 
are easier to synchronize (2). This is especially relevant for 
synchronizing clocks of sensor networks (illD . 

In Fig. [7] we plot the expander coefficient for the two mod- 
els corresponding to Fig. |6] The best parameter range is for 
small qA,B- In contrast with the distance measures, the maxi- 
mum E for Model B is network size dependent and does not 
always occur for an intermediate <7B-value. Model A, on the 
other hand, has a maximum for small g'^-values. For both 
models there is a small qa^b value where E is high and stable 
(indeed accentuated) as the system grows. This stable value is 
higher for Model A, which indicates that Model A can achieve 
higher performance than Model B for the same system size. 
For intermediate ^-values, E seems to converge for Model 
B, whereas it decreases (roughly) logarithmically for Model 
A. Because of this, in implementations where the system size 
is unknown a priori, and the model parameter is hard to con- 
trol exactly. Model B might be advantageous. Model A also 
shows a second incipient peak for the larger sizes, much lower 
than the low-^yi peak and probably uninteresting for practical 
applications. 



F. Robustness 

In addition to resilience to worst case attacks measured by 
the spectral gap, it is desirable that network graphs are robust 
to random failures. This is, in fact, similar to demand that the 
turnover of agents should not worsen our objective functions 
much. To investigate how random failures affect the network 
we measure Ad — the difference in average pathlength if a 
random node is disconnected from the network. From Fig. [8] 
we can see that this quantity is also optimized at low Qa.b- 
values. But in contrast to the behavior of E,qA -Q is optimal 
for Model A and a small qs (around - 0.2) is optimal for 
Model B (A^ = 1000). Like other quantities. Model B is less 
^g-dependent than Model A is ^^-dependent. Both models 
have similar magnitude of Aj. 
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IV. SUMMARY AND CONCLUSIONS 

We have investigated the optimization of network topol- 
ogy of the medium access layer of wireless networks, and 
found that agents with simple but heterogeneous rules for 
finding neighbors are capable to form network with shorter 
average distances, better expander properties and relatively 
high robustness (compared to random, or Unit Disc Graph ap- 
proaches). The two models we test have parameters control- 
hng the ratio of long-range interactions. We find that close to 
(but not at) one end of parameter space (where most connec- 
tions are short) is the region where optimal topologies exist. 

For most score functions — average distance, diameter and 
spectral gap — our Model A (where agents attach to spa- 
tially close others with a probability qa, otherwise to random 
agents) gives better values than our Model B (where there 
are two classes of agents, one making only short connections, 
one making only random connections). The good topological 
properties of Model A are however outweighed by the eff'ect 
of fast decreasing average degree as qa increases. Model B, 
on the other hand, has the advantage that its performance is 
less dependent on the parameter value. 
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